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Abstract 

The  Helmholtz  integral  equation  may  be  formed  by  combining  the  scalar  wave  equation 
with  Euler's  equation  for  motion  within  a  fluid.  The  solution  of  this  integral  equation 
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Sound  Wave  Scattering  from 
a  Rigid  Sphere 

1.  Introduction 

Anechoic  coatings  when  applied  to  objects  have  the  ability  to  absorb  and 
modify  the  scattering  behaviour  of  incident  sound  waves.  The  first  step  in  the 
investigation  of  anechoic  coatings  is  to  consider  their  effect  when  applied  to 
simple  objects  such  as  a  sphere  or  spherical  shell.  Before  an  analysis  of  this  is 
possible  it  is  necessary  to  have  a  firm  grasp  of  the  scattering  characteristics  of 
the  body  prior  to  the  coating  being  applied. 

The  purpose  of  this  report  is  to  investigate  the  scattering  characteristics  of  a 
rigid  sphere.  The  mathematical  derivation  of  the  scattered  pressure  using  the 
Helmholtz  integral  equation  essentially  follows  that  given  in  the  earlier  work  of 
Junger  and  Feit  (1].  When  this  is  combined  with  various  scattering  parameters 
found  in  the  literature  {2-5J,  such  as  target  strength,  reflection  factor  and 
intensity  field  around  the  body,  the  scattering  behaviour  for  the  sphere,  may  be 
characterized.  It  will  be  necessary  to  extend  this  work  later  to  consider  simple 
elastic  spheres  and  shells  as  well  as  those  having  complete  or  partial  viscoelastic 
(anechoic)  coatings  attached  to  their  outer  surface. 


2.  Helmholtz  Integral  Equation 

2.1  Introduction 

A  sound  wave  can  be  thought  of  as  a  time  dependent  pressure  fluctuation  P 
around  the  static  pressure  P,  in  a  compressible  fluid,  such  as  water.  A  sound 
field  is  generated  by  a  vibrating  elastic  structure  in  contact  with  the  fluid,  and  a 
sound  field  can  be  modified  by  the  presence  of  an  object  that  reflects  or  scatters 
the  incident  sound  waves.  The  temporal  and  spatial  variations  of  the  pressure 


fluctuation  are  governed  by  the  wave  equation 


V2P 


_L  dp2 
71F 


(2.1) 


where  c  is  the  wave  propagation  speed  through  the  fluid.  If  we  restrict 
ourselves  to  steady-state  situations  associated  with  pressures  varying 
periodically  with  time,  then  any  function  <tKf)  which  is  periodic  with  penod  T, 
can  be  represented  in  the  form  of  a  Fourier  series 


<I>(f)  *  -ia0  ♦  E  [ancos(ncof)  ♦  bnsin(nwf)l 

2  n*1  (22) 

=  E  [cn(e'ltD,)n] 

n«-~ 


where  0  =  2it  IT  is  the  fundamental  frequency.  The  coefficients  are  computed 
using  the  relations 


T  T 

an  =  1.  J  0 (t) cos (nat)dt,  bn  =  L  J  0(t)sin(nat)dt 

o  o 


(23a) 


T 

c„  =  1  j’d>(f)(eK0,)ndf  (2-3b) 

o 


In  this  report,  only  linear  systems  are  considered,  where  the  11th  term  of  the 
response  is  associated  exclusively  with  the  term  of  the  excitation. 

Therefore,  no  restriction  is  imposed  on  the  range  of  applicability  of  the  solution 
if  one  considers  a  harmonic  excitation  involving  a  single  frequency  only.  The 
time  dependence  of  the  wave  equation  ma>  be  factored  out  by  inserting  into  the 
equation  a  solution  of  the  form 


P(r,f)  =  P(  r)e-‘wt 


(2.4) 


to  yield 
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(25) 


(V2  ♦  k2)P(i)  »  0 

This  equation  is  referred  to  as  either  the  "Helmholtz  equation”  or  "steady-state 
wave  equation",  where  the  acoustic  wavenumber  is  given  by  it  =  (B/c).  Euler's 
equation  of  motion  within  a  non-viscous,  homogeneous  fluid  is  given  by- 


VP  = 


(2.6) 


where  p  and  v  are  the  density  and  velocity  respectively.  For  a  non-viscous 
fluid,  pressure  can  only  exert  forces  in  a  direction  normal  to  a  surface  boundary 
Sy  and  the  above  equation  becomes 


VP  •  n 


( v  •  n) 


on  S0 


(2.7) 


where  n  is  the  unit  outward  normal  to  S0.  Letting  the  normal  component  of 
velocity  be  given  by 


zcdr^f)  = 


where  rQ  defines  the  surface  S<y  results  in 


w 


dw 

7T 


( -  m)w 


and  the  above  boundary  condition  becomes 


dP 

7^ 


(icop)u'  on  SQ 


(2.8) 


An  elastic  body  vibrating  with  a  surface  velocity  w  will  produce  a  pressure  field 
of  the  form  and  magnitude  given  by  the  homogeneous  equation  (2.5),  and  the 
non-homogeneous  boundary  condition  (2.8).  Therefore,  equations  (2.5)  and 
(2.8)  may  be  regarded  as  defining  a  pressure  radiation  problem.  Instead  of 
solving  these  directly,  it  will  be  easier  to  combine  these  two  equations  into  one 
expression,  to  form  the  Helmholtz  integral  equation,  which  may  then  be  solved 
using  the  Green's  function  technique. 
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2.2  Theory  of  Green's  Functions 


The  Green's  function  technique  may  be  explained  as  follows.  The  equation  for 
a  pressure  field  in  the  presence  of  a  source  is  an  inhomogeneous  partial 
differential  equation  of  the  form 


(V  *  k2)P  =  q  (~9) 

where  q  is  the  source  density,  i.e.  the  strength  of  the  source  per  unit  volume. 
Conversely,  4  is  set  to  zero  in  the  above  equation  for  a  field  with  no  sources 
present.  Boundary  conditions  requiring  the  field  to  be  zero  at  the  surface  are 
termed  homogeneous  boundary  conditions.  For  non-homogeneous  boundary 
conditions,  the  boundary  values  may  be  said  to  be  caused  by  a  surface  layer  of 
sources,  corresponding  to  an  imhomogeneous  equation.  Therefore  in 
summary,  when  either  the  equation  or  boundary  conditions  are  inhomogeneous, 
sources  may  be  said  to  be  present.  When  both  are  homogeneous,  no  sources 
are  present. 

Green's  method  is  obvious  enough  physically.  To  obtain  the  field  caused  by 
a  distributed  source,  we  calculate  the  effects  of  each  part  of  the  source  and  add 
them  all.  If  G(r  |  r0)  is  the  field  at  the  observer's  point  r  caused  by  a  unit  point 
source  at  r0,  then  the  field  at  r  caused  by  a  source  distribution  q(r0;  is  the 
volume  integral  of  G(r  |  r0)^(r0)  over  the  whole  range  of  r0  occupied  by  the 
source.  G  is  called  Green's  function  and  is  a  solution  of  the  inhomogeneous 
Helmholtz  equation 


(V2  +  k2)G( r|rQ)  =  6(r 


(2.10) 


where  6^-^,)  is  the  Dirac  delta  function.  G^rl^)  satisfies  the  reciprocity 
principle 


cedi) 


(2.11) 


Boundary  conditions  can  be  satisfied  in  the  same  way.  We  compute  the  field 
at  r  for  the  boundary  value  P0  (or  normal  gradient  dP/dna )  as  zero  at  every 
point  on  the  surface  except  for^,s  (which  is  on  the  surface).  At the 
boundary  value  has  a  delta  function  behaviour,  so  that  its  integral  over  a  small 
surface  area  near  r0s  is  unity.  This  field  at  r  (not  on  the  boundary)  we  can  call 
G(t,LEoS)P0  (°r  g(lI.5jS)  3P/3n0)  over  the  boundary.  G  here  is  also  Green's 
function. 

One  can  solve  the  "inhomogeneous  equation"  for  a  field  caused  by  a  source 
distribution  q,  by  means  of  a  product  of  the  source  density  with  a  Green's 
function  integrated  over  a  volume.  The  solution  of  a  homogeneous  equation 
having  specified  values  on  a  surface  can  be  obtained  in  terms  of  another 


Green's  function,  integrated  over  the  boundary  surface.  This  Green's  function 
happens  to  be  identical  to  the  first.  Therefore,  Green's  function  for  a  radiation 
problem  governed  by  the  equations 


(V2  +  k2)P  =  0,  ££  =  (tup) w  on  S 

Bn 


(2.12) 


may  be  determined  from 


(V2  *  kz)G(. rjr,)  -  S(r -r^,  !gL  =  0  onS0  (2.13) 


Before  beginning  to  formulate  the  Helmholtz  integral  equation,  there  remains 
the  task  of  investigating  the  behaviour  of  Green's  function  G(r  |  r0)  for  the 
situation  where  observation  point  r  and  the  source  point  r^  are  such  that  the 
magnitude  R  -  i  r  -  r0 1  is  small  compared  with  the  distance  of  either  point  from 
a  surface  boundary.  In  this  case,  G  is  just  a  function  of  R,  i.e.  the  source  is 
completely  symmetrical,  so  that  G  cannot  depend  on  the  direction  of  R,  only  on 
its  magnitude.  This  Green's  function  is  termed  the  free  space  Green's  function 
g,  and  has  a  singularity  at  R  -  0.  Therefore,  for  a  general  problem  involving  a 
surface  boundary,  G(r  |  r0)  may  be  separated  into  two  parts:  a  part  F  which 
depends  on  the  boundary  conditions  at  and  a  part  g,  which  is  continuous 
everywhere  except  at  r  =  rQ  and  which  is  a  function  of  R  alone.  In  the  next 
section,  the  free  space  Green's  function  for  a  three  dimensional  situation  is 
derived. 


2.3  The  Free  Space  Green's  Function 


To  begin  with,  the  free  space  Green's  function  must  be  a  solution  to  the 
inhomogeneous  Helmholtz  equation 


(V  ♦*2)g<|r-.rD|>  »  Sfr-r.) 


(2.14) 


To  construct  a  solution,  a  spherical  system  of  coordinates  is  chosen  where  the 
origin  coincides  with  the  singularity  in  equation  (2.14).  Thus  the  distance 
|jr  |  reduces  to  the  radial  coordinate  R.  In  this  spherical  system,  the 
solution  is  a  function  of  R  only  and  not  of  the  angular  coordinates.  Therefore, 
except  at  the  origin.  Green's  function  satisfies  the  homogeneous  Helmholtz 
equation 
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a2  +  23 

IF  ’EM 


g(R) 


0,  R  >  0 


(2.15) 


which  has  the  general  solution 


g(R)  =  (l/RHA:^  +  Be _ikR> 


(2.16) 


Earlier  in  this  report,  the  Helmholtz  equation  was  obtained  from  the  wave 
equation  by  specifying  the  time  dependence  as  simple  harmonic,  i.e. 

Pf&t)  =  P[ r)e'i0*.  If  g  is  to  be  used  as  a  solution  of  the  wave  equation,  then  the 
complete  solution  is 


g(R)e-iat  «  (l/j?)[V(kR'“,>  ♦  Be*l0‘R*wt)] 


To  obtain  waves  travelling  outward  from  the  source  point,  we  must  set  B  to 
zero,  such  that 


g{R)e'i(at  =  (A/R)e,(kR'“t) 


represents  an  outgoing  wave.  The  coefficient  A  may  be  determined  by 
integrating  equation  (2.14)  over  a  spherical  volume  element  of  radius  a 
concentric  with  the  singular  point  R  =  0,  and  then  taking  the  limit  as  a  -»  0. 
This  process  yields 


lim 

a-*0 


SS  $v2*iv 


1 


From  Gauss'  divergence  theorem,  this  may  be  rewritten  as 


lim 


-  "I 


where  n  is  a  unit  vector  that  points  into  the  volume  V.  Now 


-n 


-dg 

IT 


A 

^  as  a  — >  0,  R  *  a 

a" 


Since  this  formulation  is  independent  of  the  angular  coordinates,  the  surface 
integral  is  obtained  simply  by  multiplying  this  result  by  the  area  of  the  sphere, 
W ,  which  yields  A  ~  -  l/4rt.  Substituting  for  A  and  replacing  the  spherical 
coordinate  R  by  |,r  - i ,  the  free  space  Green's  function  is  obtained 


-e 


(2.175 


2.4  Formulation  of  the  Integral  Equation 

In  order  to  formulate  the  Helmholtz  integral  equation,  we  will  require  the  use 
of  Green's  theorem  which  may  be  derived  as  follows.  From  Gauss'  divergence 
theorem 


n  into  the  volume  V  (2.18) 

Consider  two  scalar  functions  U  and  V  and  let  the  function  F  be  given  by 

F  -  UVV-VVU  (2.19) 


Substituting  for  F  in  the  above  integral  expression  yields  Green's  relation 
between  volume  and  surface  integrals,  i.e. 


*  -JJlUVV- VVU1  •  dS  (22m 


To  formulate  the  Helmholtz  integral  equation,  consider  the  two  inhomogeneous 
partial  differential  equations 
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<V*  ♦*:*)?( r)  =  ‘(r) 


(2.21) 


(V2**2^^) 


S(£“-£o> 


(222) 


Multiplying  equation  (2.22)  by  P( t)  and  equation  (221)  by  G'rJ^  and 
subtracting,  exchanging  Horn,  at  the  same  time,  yields 


Integrating  over  the  volume  of  S^,  and  using  the  property  of  the  delta  function, 
namely 


[  Pit) 


jr  inside  V 
i  outside  V 


(223) 


yields 

p(r)  *  / J/l (p^)V2GCr!^)  -G(r|rJV2P^)]dU  ./J/fGftl^)^)]^ 

(224) 

The  first  term  may  be  expressed  as  a  surface  integral  by  using  equation  (220) 
with  the  substitutions  U  =  PtrJ,  V  =  G^r  Ij^),  to  yield 

pW  =  -J/fp^>VG'tl£o)-G(r|^)VP^))-dS*/J/{G(ri^)q(r0)lc'V 

(225) 

Earlier  in  the  report,  we  had  the  equations  for  the  radiation  problem  as  being 

(V2*fc2)P  ^  0 ,  ££  =  (i'(op)w  on  SQ,  q  =  0  (226) 

on 


Substituting  equation  (2.26)  into  (2.25),  yields 
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P(r)  ~  -Jj  [Pfr,,)  -5—  " (»’<0p)G(£|ro)a>(jo)]rfS,  i  =  S0,  r  inside  V 


We  have  now  come  up  against  a  problem,  the  solution  to  the  above  expression 
represents  the  pressure  fluctuation  Pii )  inside  a  boundary  surface  Sv  yet  we  are 
trying  to  formulate  an  equation  for  the  radiated  pressure  outside  a  vibrating 
body.  To  overcome  this  problem,  the  volume  V  inside  the  boundary  S  is 
defined  as  the  volume  V  between  the  vibrating  surface  S0  and  that  of  an  infinite 
sphere  of  radius  r,.  The  expression  for  the  pressure  field  then  becomes 


^  -<J‘0P>G(r;Lro)i4;(ro))<iS,  _r  inside  S  =  S0  ■ 


where  Sj  is  the  surfr  e  of  the  infinite  sphere.  Unless  the  surface  integral  at 
infinity  in  the  above  equation  vanishes,  one  reaches  the  conclusion  that  the 
pressure  field  is  not  uniquely  determined  by  the  boundary  conditions  over  the 
radiating  surface  S0.  The  surface  integral  over  S1  may  be  written  as 


r.  ,  .  dP , 

-  G(r|r. )  _) 
~  dn 


um  2, ,  ,  dG  r,  ,  .  df. 


Because  there  are  no  boundary  conditions  on  the  infinite  sphere,  the  Green's 
function  G  may  be  replaced  by  the  free  space  Green's  function,  i.e. 


eflr-r-l)  **  _ 

3  ~  An  i 


The  surface  integral  becomes 


Um  r .,  D/  .  3P, 

rj— ko  ri  ItkPhj)  - 


To  satisfy  this  condition,  the  pressure  must  decrease  with  increasing  r1  ns 
If]  -  r0 1  **  or  faster.  According  to  Jungc -  and  Feit  [1],  this  is  the  case  for  all 


sources  of  finite  extent,  such  as  the  sphere.  Therefore,  the  integral  equation 
can  be  wntten  as 


PC r)  -  -  -  <ia>p)G^!^)ix'(r0)JiiS(r0),  £  outside  Sa 


(228) 

which  is  typically  known  as  the  Helmholtz  integral  equation.  This  integral 
may  be  reduced  to  a  simple  expression  involving  only  Green's  function  and  the 
known  surface  velocity  on  the  boundary,  providing  that  Green's  function  can  be 
constructed  in  a  way  that  satisfies  the  Neumann  boundary  condition 


8G 

IK 


0  on  the  boundary  Se 


(229) 


then  P(r)  reduces  to 


P(r)  =  (/tap)  j  J"  f  GCr^rr^)  ti?  I  dS0 


(2.30) 


i.e.  the  field  is  just  the  surface  integral  of  Green's  function  multiplied  by  the 
boundary  values  as  mentioned  in  section  2.2.  If  instead  of  a  surface  layer  of 
sources,  a  distributed  source  4  is  the  only  source  present,  then  from  equation 
(225),  the  pressure  field  reduces  to  the  volume  integral 


P(r)  =  JJJ[G(r|rQ)^:  )\dV 


(231) 


as  mentioned  previously,  provided  that  G  satisfies  the  Neumann  condition  on 
the  boundary  for  this  case  also. 


3,  Pressure  Field  for  a  Spherical  Radiator 

3.1  Green's  Function  for  a  Vibrating  Sphere 

In  order  to  obtain  a  solution  to  the  Helmholtz  integral  equation  for  a  vibrating 
source.  Green's  function  G  must  be  obtained.  As  mentioned  earlier  in  section 
22,  a  general  Green's  function  may  be  constructed  by  adding  a  term  T  to  the 
free  space  Green's  function,  where  G  satisfies  the  boundary  condition 
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(3.1) 


dC 

~5n 


0  on  Sg 


This  additional  term  contains  the  natural  eigenfunctions  for  a  sphere  (see 
Appendix  A  for  eigenfunction  derivation).  The  first  step  in  the  construction  of 
Green's  function  is  to  express  the  free  space  Green's  function  in  spherical 
coordinates.  To  do  this,  we  must  first  separate  g  as 


*  - 


.  -it 

4it|_r-r0(  4  it 


costfc  Ijl-JEoI)  tsm{k 


*l£T£oi 


(32) 


From  trigonometry,  the  distance  between  the  field  point  r  and  the  source  point 


Ijr-iol  =  (r2  +rl  -2rr0cosv)1/2 


(33) 


where  y  is  the  angle  between  r  and^.  The  components  of  g  can  be 
represented  in  terms  of  y  by  the  expressions  [6] 


c°s(*l£-fJ)/*l£-r»l  *  ~  £  yn(fcr)(2«*l)/n(fcr0)Ptl(cosMr) 


(3.4) 


I£t£oI)  /  *  l£-£ 


-JO' 


E  ;n(W(2n  +  l);n(fcr0)Pn<cosV) 
n^o 


where  jn(kr),  yn(tr0)  are  spherical  Bessel  functions  and  Pn(cos  y)  is  the  Legendre 
polynomial.  Thus 


£<U,-.£oP  *  ~(ik/4n)  £  (2n  *l)jTi(kr0)hn(kr)PTx(cosy) 

n»o 


(33) 


where 


ha(kr)  =  /„<*»•) +«y n<*»0 


(3.6) 
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Due  to  symmetry  about  the  d  =  0  axis,  the  angle  y  may  be  written  as 
y  =  t)  -  d0<  and 


Pn(cosy)  =  Pn(cosd)Pn(costi0) 


Thus,  the  free  space  Green's  function  in  terms  of  (r,  0)  and  (r^  i30)  instead  of  r 
and  r„  is 


Sdi-rJ)  -  -(ik/4n)  £  (2n ♦l)Pn(cosri)Pn(cosd0);n(^0)hn{fcr)  <3.7) 

ii«o 


The  complete  Green's  function  may  be  represented  as 

G(£(rp)  -  X<l£-^£oP  +  nX!io)  a8) 


and  must  satisfy  the  boundary  condition 


ar  dg 


at  rQ  -  a  (Sphere  surface) 


(3.9) 


Let  the  (r,  r0)  modal  components  of  g  and  T  be  represented  by  g„  and  Fn 
respectively.  Assuming  that  Tn  is  made  up  of  the  eigenfunctions  of  the  sphere 
and  must  also  satisfy  the  redprodfy  condition  Tn(r  |  r0)  =  rn  (r0  |  r)  we  then  can 
expect  f„  to  have  the  form 


fn(r,r0)  *  Ahn{kr0)hn(kr) 


(3.10) 


where  from  equations  (3.9)  and  (3.10) 


9rn/dr0  =  kAh^ikr  0)hn(kr)  =  -kj^{ka)hn(kr) 


Thus 


f'n<r'ro>  *  - hn (kr)hn ikra ); '  (ka)/h!tt  (to) 


(3.11) 
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Adding  gn  and  fn  and  using  the  relation 


fn(ka)y^(ka)  -yn  (&)/„  (ka)  *  (kaf2 


(3.12) 


yields 


in  *  Fn 


(kafh^Oa) 


at  rQ  *  a 


(3.13) 


The  function  T  has  the  same  (d,  dependence  as  g.  Therefore,  the  resulting 
Green's  function  for  a  spherical  source  may  be  written  as 


C(r,tf!a,a0) 


r 


(2  n  +DPn(cosd)Pn(cost30) 


hn(kr)  i 


(3.14) 


3.2  Solution  to  Integral  Equation 

Now  that  Green's  function  has  been  determined,  the  radiated  pressure  field 
may  be  obtained  from  the  Helmholtz  integral  equation  (2.30),  ie. 

P(r)  =  (tap)  fflGirlr^uHr^JdS 

~  (3.15) 

=  (iap)(2jrtf2)  J [G(r,0  |fl,i30)a/W0)]sini30dd0 


To  solve  this  integral,  it  is  helpful  to  express  the  surface  velocity  as 


ut to)  =  lVnPn(cosi3) 
n»o 


(3.16) 


where  the  coefficients  of  the  series  are  given  by 


cosd 


(3.17) 


1 

%  =  —  *!-  J  Pn(TDw(T))iT1,  T1  = 
2  -1 


Using  the  relation 

l 

/  =  (26nJn)/(2n>l) 

-l 


results  in  the  final  expression  for  the  pressure  field,  ie. 


P(r/d)  *  i(pc)  £  WnPn(cosd) 

n«o 


h^ikr) 

7 


(3.18) 


(3.19) 


4.  Scattering  of  a  Plane  Wave  from 
a  Rigid  Sphere 

Consider  a  distant  point  sound  source,  which  generates  a  continuous  sound 
wave.  Far  away  from  this  source  and  over  suitably  restricted  regions,  these 
waves  may  be  said  to  approximate  plane  waves.  Consider  these  plane  waves 
to  be  incident  upon  a  stationary  rigid  sphere,  whose  surface  SQ  is  given  by  the 
position  vector  shown  in  Figure  1.  If  we  define  a  wavenumber  vector  k, 
whose  magnitude  is  k  and  which  lies  in  the  direction  of  wave  propagation,  then 
the  incident  pressure  wave  may  be  written  as 


P[(£,f)  =  P0expi(k-r -cat) 

In  spherical  coordinates,  a  wave  incident  from  the  d  =  180°  direction  is 
represented  by 

Pj(r,d)  *  P0expi(fcrcosd)  (4.1) 

The  introduction  of  the  rigid  boundary  Sq,  produces  a  disturbance  in  the 
pressure  field  P,,  which  is  termed  the  scattered  pressure  where  the  <~ 
symbol  refers  to  the  infinite  impedance  of  the  scattering  surface.  The  resultant 
pressure  in  the  presence  of  the  scatterer  is  defined  as 


Pit)  -  Pt(r)  ♦  P5Jr> 


(4  2) 


Scattered 
sofencal  wave 


Figure  1:  Sound  scattering  from  a  rigid  sphert. 

In  order  to  calculate  the  scattered  pressure  and  hence  the  resultant  pressure 
field,  we  need  to  make  use  of  the  idea  of  re-radiation.  This  amounts  to 
regarding  the  boundary  surface  as  having  a  distribution  of  virtual  sources, 
whose  nature,  strengths  and  phases  are  determined  by  the  incident  wave  and 
by  the  properties  of  the  water  and  the  boundary  medium.  The  radiation  from 
these  virtual  sources  then  determines  the  behaviour  of  the  reflected  or  scattered 
wave.  To  use  this  method,  we  need  to  investigate  the  components  that  make 
up  the  velocity  distribution  on  the  surface  of  the  sphere.  Since  the  boundary  of 
the  sphere  is  rigid,  the  resultant  velocity  on  S0  normal  to  the  surface  must  be 
zero,  i.e. 

wit)  =  Wj(r)  +  w^tr)  *  0 

The  resultant  surface  velocity  may  be  related  to  the  pressure  gradient  via 
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-j—  =  (icop)w  on  SQ 


(4.4) 


which  may  be  rewritten  as 


•to, 


J_  9Pi 
cap  "SF 


_ L  (k  -n )P-.  on  SD 

cap 


(45) 


Therefore,  may  be  regarded  as  the  radiated  pressure  field  due  to  a  vibrating 
spherical  source  with  the  velocity  distribution  and  hence  the  solution  to 
the  Helmholt2  integral  equation  for  a  spherical  radiator  may  be  used.  In  order 
to  apply  the  Helmholtz  solution,  the  velocity  must  first  be  expressed  in 
terms  of  Legendre  polynomials  as 


*  E  WnPn(cosd) 

n-o 


To  do  this,  we  proceed  as  follows.  The  expression  for  a  plane  wave 
Pj<r,d)  =  PQexpi(kra  d) 


can  be  written  in  terms  of  a  series  of  concentric  spherical  waves  as 


P^r.-d)  -  P0  £  (2n  ♦  1  )(DnPn  (cos  d)/n  (lx) 
n*o 


(4.7) 


The  coefficients  for  then  become 


Wn  «  <iP0/pc)(2n+l)(i)nj'(*a) 


(4.8) 


From  the  solution  for  a  spherical  radiator,  the  radiated  pressure  field  is 
represented  by 
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P(r,d)  *  (ipc)  VVnPn(cosi3)/in  (kr)  /h^  (ka) 
n«o 


Substituting  equation  (4.8)  for  Wn  into  equation  (4.9),  yields  the  rigid  scattered 
pressure  field 


P^(rid)  *  -P0  £  <2»«-l)(j)nPn(cost»/'  (ka)hn(kr)/h„  (ka)  <4-10) 

n«o 


The  above  series  may  be  simplified  somewhat  if  we  restrict  ourselves  to  large 
distances  r  from  the  boundary  such  that  kr  »  1.  Then  the  Hankel  function 
k^kr)  may  be  approximated  by 

(l/fcr)eikr(i)'(n'>1)  kr>>n2*  1  t4>11) 

Substituting  h^Ocr)  into  P^  yields  what  is  termed  in  this  report  as  the  "far-field" 
scattered  pressure 


P^frifi)  -  dP0/kr)exp(ikr)  £  (2n  +l)Pn(cosfi)/'  <*»>/&'  (ka)  (412) 

n»o 


"Near-field"  in  this  report,  means  that  equation  (4.10)  is  used  for  the  rigid 
scattered  pressure,  while  for  "far-field"  (i.e.  kr  »  1),  equation  (4.12)  is  applied. 
The  first  ten  terms  of  the  series  have  been  found  to  be  sufficient  to  determine 
the  scattered  pressure  in  the  ka  range  under  investigation. 


5.  Methods  of  Presenting  the  Scattering 
Solutions 


In  this  chapter  mathematical  expressions  are  developed  for  a  number  of 
parameters  which  are  useful  in  examining  the  scattering  behaviour  of  the  rigid 
sphere.  These  parameters  include  (i)  the  resultant  pressure  field  magnitude  in 
the  backscatter  direction,  (ii)  the  target  strength  of  the  scattered  (iii)  the 
reflection  factor  around  the  scattering  body  and  (iv)  the  intensity  field  around 
the  scatterer.  To  determine  these  various  parameters,  the  following  values  for 
the  propagation  medium  will  be  used. 


23 


Water 


Density  of  medium  p  =  1000.0  (kg/m3) 

Speed  of  sound  c  -  1500.0  (m/s) 


5.1  Resultant  Pressure  Field  in  Front  of  Scatterer 

Here,  the  resultant  pressure  magnitude  as  a  function  of  distance  from  the 
scatterer  is  computed,  in  the  direction  towards  the  incident  waves.  The 
normalized  resultant  pressure  field  for  the  rigid  sphere  may  be  written  as 


P  =  |Pj/P0  -  KPi/Po)  MP^/P^I 


(5.1) 


where  the  incident  and  rigid  scattered  pressures  dose  to  the  sphere  have  been 
determined  previously  to  be 


and 


Pj/P0  *  P,N  (r,d)  /P0  =  exp(iftcosd) 


(5.2) 


p£>o  =  P^r.W/Po 


oo 

-r 


(2n  +l)(i)nPn(cosO) 


hn(ir) 


(53) 


where 


x  -  ka,  f  -  r/a 


(5.4) 


5.2  Target  Strength 


Target  strength  is  a  far-field  parameter  which  for  convenience  is  defined  as  the 
ratio  of  the  scattered  intensity  considered  at  one  metre  to  the  incident  intensity 
expressed  in  dedbels,  i.e. 


/_(r,d) 

TS  =  10  logj0  -22. - jr„ 

10 


(55) 


24 


where  f0  is  the  intensity  of  the  incident  (plane)  wave  as  defined  in  equation 
(5.13).  The  scattered  intensity  is  proportional  to  the  square  of  the  pressure,  and 
in  the  far-field  is  proportional  to  1/r.  Therefore,  the  target  strength  may 
rewritten  in  the  form 


TS  -  20  logjg 


rPjL(r,-d) 

\-Jl _ 

P 

o 


(5.6) 


where  PFSOQ  is  the  scattered  pressure  in  the  far-field  and 


(r/P0)P^(r,i3)  -  i(a /x)exp(izr)  £  (2«  +l)Pn(cosi3)[ /„  (t)J 

n-o 


Results  obtained  for  different  fca  values,  will  be  compared  with  the  commonly 
used  expression  for  the  target  strength  of  a  sphere,  i.e. 


TS  =  10  log10  (a2/4) 


(5.8) 


5.3  Reflection  Factor 


Using  the  expressions  in  the  previous  section,  the  ratio  of  scattered  to  incident 
pressure  may  be  written  as 


f 

a 

' 

gitf(l-eo*e) 

7> 

(:/x)  £  (2n+l)Pn(cosfl) 


/n<*> 


KM 


(5.9) 


The  1/r  term  characterizes  the  spherical  spreading  of  the  scattered  wave,  whilst 
the  exponential  term  takes  into  account  the  phases  of  the  incident  and  scattered 
waves.  The  reflection  factor  can  be  made  independent  of  these  terms  by 
defining  it  as 


R  =  (r  /a)  j  e  ~ixKl  ■cos*>  [  p£,  /Pj  ]  | 
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to  yield 


R  = 


|  (»7t)  5^  (2n  +  l)Pn(costi) 
11*0 


/> 

^(T) 


(5.10) 


The  reflection  factor  will  be  computed  as  a  function  of  direction  around  the 
body,  for  various  ka  values. 


5.4  Intensity  Field  Around  Scatterer 

Using  the  solution  to  the  resultant  pressure  field,  the  construction  of  two 
different  types  of  time-averaged  intensity  plots  is  possible. 

(A)  Plots  of  the  intensity  magnitude  around  the  body 
(both  near-field  and  far-field) 

(B)  Plots  of  the  intensity  vector  field  around  the  body 
(near-field  only) 

For  both  types  A  and  type  B  plots,  the  intensities  for  the  scattered  and  resultant 
pressure  fields  have  been  computed.  The  radial  and  transverse  components  of 
the  intensity  vector  (averaged  over  one  period  of  the  wave)  are  defined  as 


IAr,-Q)  =  1  fk[P(r,t5)iir*<r,fl)J  (5.11a) 

r  2  r 

=  1.  RetP(r,d)u#(r,iJ)]  (5.11b) 


where  Re  denotes  the  real  part,  *  the  complex  conjugate  and  Uj,  u#  the  velocity 
components  corresponding  to  the  pressure  field  P,  whether  it  be  a  scattered  or 
resultant  pressure.  The  magnitude  /  and  direction  p  of  the  intensity  vector  are 
given  by  the  formulae 


I  =  l420'.»W.2('.<01ly3 

♦  =  0  *arctan[/#(r,0)//r(r,ft)] 


where 
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d  =  arctan{z/x),  r  =  (x2+z2j1/'2 

and  x,  z  are  shown  in  Figure  1.  For  the  polar  intensity  plots  of  type  A,  the 
intensity  value  is  first  divided  by  the  incident  intensity 


(5,13) 


and  then  normalized  with  respect  to  the  maximum  value  that  it  obtains  over  the 
range  of  d  varying  from  0®  to  360°.  For  the  type  B  intensity  vector  plots,  each 
vector  has  first  been  divided  by  I#  and  then  normalized  with  respect  to  the 
magnitude  of  the  largest  vector  in  the  plot.  The  square  root  of  the  magnitude 
is  then  calculated  and  multiplied  by  the  grid  spacing.  The  vectors  are  plotted 
over  the  coordinate  range  -  2  S  x/a  <  2  and  0  S  z/a  <  2,  where  a  is  the  radius  of 
the  spherical  scatterer. 

in  order  to  compute  the  intensity,  the  velocity  field  associated  with  the 
pressure  field  must  be  found.  This  may  be  determined  from  Euler's  equation 
VP  =  (l(op  hi,  which  in  component  form  becomes 


f  \ 

1 

f  \ 

-x 

dP 

(  \ 

l  i 

r-.i 

a 

(tJ 

17'  “*  ‘ 

Tc) 

,*rJ 

(5.14) 


Using  the  near-field  and  far-field  pressure  expressions  given  in  sections  5.1  and 
52,  the  following  velocity  components  may  be  obtained. 


Velocity  due  to  incident  pressure  P{ 

(l pc/P0)  [wrJt  =  cosdetocos* 

(pc/PpHii^lj  =  -sini3eilacost} 


(5.15a) 


(5.15b) 
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Velocity  due  to  rigid  scattered  pressure  P^ 
Near-field: 


(pc/P0)[iirN]^. 


=  £  (2 n  *  PUT'1  Pn(cosf»  !liX).  h‘_  (tf) 
r-o  ‘  Jl^(x) 


(5.16a) 


(pc/P0)tu^ 


“  ,  ,  j„  (x) 

=  -(sinti/xf)  £  (2n+l)(On'1Pn(cos0)  J! — An(xf) 

n-o  h ^  (T) 


(5.16b) 


Far-field: 


(pc/P^l^  --  (i/xr)enf  £  (2n+l>Pn(cosd)_l 


(5.17a) 


/!„(!) 


(pc/P^lu^J,.  =  (-sind/fcf)2?**'  £  CwDP'teosd)^!.  (5.17b) 


Vx> 


Therefore,  by  inserting  the  pressure  and  velocity  formulae  given  respectively  by 
equations  (5-2),  (53),  (5.7)  and  (5.15),  (5.16),  (5.17)  into  equation  (5.11),  the 
scattered  and  resultant  intensity  expressions  may  be  obtained. 


6.  Results  and  Discussion 

6.1  Resultant  Pressure  Field 

Figures  2(a)  to  2(d)  show  the  normalized  resultant  pressure  as  a  function  of  r/a 
for  the  ka  values  1,  2,  4  and  8.  In  these  figures,  the  reflected  wave  combines 
with  the  incident  wave  to  form  a  standing  wave  whose  scattered  component 
decreases  proportionally  to  1/r  where  r  is  the  distance  from  the  scatterer.  As 
the  frequency  increases,  (ka  increasing),  the  pressure  at  the  surface  approaches 
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twice  the  incident  pressure  amplitude,  which  is  the  case  for  a  plane  wave 
incident  upon  a  rigid  wall. 


6.2  Target  Strength 


Figure  3  shows  the  far-field  target  strength  for  the  rigid  sphere  over  the  ka 
range  from  0  to  8.  From  this  figure,  it  is  evident  that  for  frequencies 
corresponding  to  ka  greater  than  2.0  and  within  the  range  considered,  the  target 
strength  is  approximately  constant.  These  results  compare  favourably  with  the 
usual  expression  for  the  target  strength  (due  to  specular  reflection)  of  a  rigid 
sphere,  i.e. 

TS  =  10  logl0  (a2/4)  =  -6  (dB)  for  sphere  of  radius  a  -  1 


which  is  shown  by  the  dashed  line  in  Figure  3. 


r/a 

Figure  2(a):  Resultant  Pressure  in  front  of  Rigid  Sphere. 
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Figure  2(b):  Resultant  Pressure  in  front  of  Rigid  Sphere. 


Figure  2(c):  Resultant  Pressure  in  front  of  Rigid  Sphere. 
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6.3  Reflection  Factor 


Figures  4(a)  to  4(d)  show  the  bistatic  reflection  factor  around  the  scatterer  for 
the  ka  values  1, 2,  4  and  8.  From  these  figures,  we  can  see  that  for  small  ka 
values  (where  the  sphere  is  small  compared  with  the  wavelength),  the  wave 
propagates  in  all  directions.  For  large  ka  values,  part  of  the  scattered  wave 
spreads  out  more  or  less  uniformly  in  all  directions  from  the  scatterer,  and  the 
rest  is  concentrated  in  the  0°  direction.  When  the  resultant  field  is  computed 
by  adding  the  incident  and  scattered  components,  the  0°  beam  of  the  scattered 
component  interferes  destructively  with  the  unchanged  plane  wave  behind  the 
sphere,  creating  a  sharp  edged  shadow. 


6.4  Intensity 


Figures  5  and  6  show  the  scattered  intensity  magnitude  around  the  body  in  the 
near  and  far-fields  for  ka  values  1,  2,  4  and  8.  These  figures  exhibit  essentially 
the  same  behaviour  as  observed  with  the  reflection  factor.  Figures  7(a)  to  7(d) 
show  the  resultant  intensity  close  to  the  body.  The  large  side-lobes  evident  in 
Figures  7(a)  and  7(d)  arise  from  the  constructive  interference  of  the  incident  and 
scattered  waves. 

Figures  8(a)  to  8(d)  show  the  scattered  intensity  vector  field  around  the  body, 
and  illustrates  more  dearly  the  polar  intensity  plots  in  Figure  5.  Figures  9(a)  to 
9(d)  represent  the  resultant  intensity  field  around  the  body  and  show  clearly  the 
reduction  in  intensity  behind  the  sphere  as  the  frequency  increases. 


7.  Conclusion 


An  analysis  of  the  plane  wave  scattering  from  a  rigid  sphere  of  radius  a 
submerged  in  water  has  been  carried  out  and  results  are  presented  for  ka  values 
up  to  8.  Various  scattering  parameters  such  as  target  strength,  reflection  factor 
and  the  intensity  field  around  the  body  have  been  computed  to  illustrate  the 
solutions.  Resultant  intensity  vector  plots  show  dearly  the  reduction  in 
intensity  behind  the  sphere  as  the  frequency  increases. 
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Figure  4(c):  Reflection  Factor  for  a  Rigid  Sphere. 


ka  ■  8  90’ 


Figure  4(d):  Reflection  Factor  for  a  Rigid  Sphere. 
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Figure  5(b):  Scattered  Intensity  Ijlmax  around  a  Rigid  Sphere. 


Figure  5(c):  Scattered  Intensity  I/Imax  around  a  Rigid  Sphere. 


90* 


270" 

Figure  5(d):  Scattered  Intensity  Ijlmax  around  a  Rigid  Sphere. 
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90* 


r/a  -  50.0 


Figure  6(a):  Scattered  Intensity  I/Imax  around  a  Rigid  Sphere. 


270* 


Figure  6(b):  Scattered  Intensity  I/Imax  around  a  Rigid  Sphere. 
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Figure  6(d):  Scattered  Intensity  I/lmax  around  a  Rigid  Sphere. 


Figure  7(a) :  Resultant  Intensity  1/lmax  around  a  Rigid  Sphere. 


270' 

Figure  7(b):  Resultant  Intensity  1/lmax  around  a  Rigid  Sphere. 


770° 


Figure  7(d):  Resultant  Intensity  Ijlrnax  around  a  Rigid  Sphere. 
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Figure  8(c):  Scattered  Intensity  (l/ImaxP<2  around  a  Rigid  Sphere. 
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Figure  8(d):  Scattered  Intensity  [ I/Imax) 2/2  around  a  Rigid  Sphere. 
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Figure  9(a):  Resultant  Intensity  ll/lmax]lf2  around  a  Rigid  Sphere. 


Figure  9(c):  Resultant  Intensity  [Ijlmaxlir2  around  a  Rigid  Sphere. 
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Figure  9(d):  Resultant  Intensity  U/lmaxf'1  around  a  Rigid  Sphere. 
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9.  Glossary  of  Symbols 


p 

Po 

Pi 

Ps- 

w 

W1 

Ws- 

'o 

I*- 

I 

Jrlo 

c(Llic) 

gdzsob 

r,  P,  <j> 
a 

So 

n 

k 

ka 

a 

c 

P 

R 


Resultani  pressure 

Incident  pressure  wave  amplitude 

Incident  pressure 

Scattered  pressure  from  a  rigid  body  (infinite  impedance) 
Normal  velocity  component  of  resultant  pressure  at  S0 
Normal  velocity  component  of  incident  pressure  at  S0 
Normal  velocity  component  of  rigid  scattered  pressure  at  SD 
Average  intensity  magnitude  of  the  incident  pressure 
Average  intensity  magnitude  of  the  scattered  pressure  (rigid 
body) 

Average  intensity  magnitude  of  the  resultant  pressure 

Position  vectors  of  the  field  and  source  points 

Green's  function 

Free  space  Green's  function 

Spherical  coordinates 

Radius  of  sphere 

Boundary  of  sphere  (r  =  a) 

Unit  outward  normal  to  boundary  Sa 
Acoustic  wavenumber 
Dimensionless  frequency  parameter 
Angular  frequency 
Speed  of  sound  in  water 
Density  of  water 
Reflection  factor 


46 


Appendix  A 

Eigenfunctions  of  a  Sphere 


In  this  section  the  eigenfunctions  of  a  sphere,  which  are  required  in  order  to 
construct  a  suitable  Green's  function  for  a  spherical  radiator,  are  derived.  The 
Helmholtz  equation 


(V2+fc2)P  *  0 


(A.l) 


may  be  expressed  in  spherical  coordinates  (r,  0,  <h)  as 


1  d 

7-57 


(  \ 

2  ap 


1 


a 

sind 


sind 


3P 

"37 


1 


k2P 


0 


(A  2) 


where  it  is  assumed  that  the  sphere  in  question  is  symmetrical  about  the  d  =  0 
axis,  such  that  the  $  dependence  may  be  omitted.  Letting  the  pressure  be 
represented  by  P(r,  d)  =  R(r)T( d),  modifies  the  Helmholtz  equation  to 


1  d 

R  Hr 


.2  dR 
IF 


*  k2r2 


-1  1  d 

T  sind  dd 


fsind 


(A3) 


Each  side  of  this  equation  must  equal  the  same  constant  C,  i.e. 


1  d 
?Tr 


(  \ 

2  dR 

/ 

l  ir J 

\ 

0 


(A.4) 


(  \ 


1  d 

sind  7d 


J^sind 


iT 

3d 


*  CT 


0 


(A3) 


Letting  q  =  cos  0,  C  =  n(n  +  1),  Legendre's  differential  equation  is  obtained 
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(A.6) 


(i-n2)i!l 

dr? 


2t\£L  *  n(n  +1)7  =  0 

i  T) 


which  has  the  Legendre  polynomial  solution 


Tn<»>  =  Pn(cosi3) 


(A  ,7) 


The  equation  for  the  radial  factor  R  can  be  modified  by  substituting  z  =  for,  to 
yield  the  equation  for  the  spherical  Bessel  function 

z2ilj  *2z^L  +[z2-n(n  +  l)lK  *  0  <A.8) 

dz2  « 


from  which  the  possible  solutions  are  jn(fcr),  yn<lr),  nd  fcj  (2)(*cr>,  where 

/n(Jtr)  and  yn(ir>  are  the  spherical  Bessel  functions,  and  hn<l>()cr),  hj2)(kr)  are  the 
Hankel  functions  of  the  first  and  second  kind  respectively-  The  chosen 
solution  is 


R(r)  -  h"\kr) 


(A.9) 


due  to  the  fact  that  the  spherical  Hankel  function  hna)(kr)  corresponds  to  an 
outgoing  wave,  appropriate  for  the  situation  where  acoustic  energy  is  being 
radiated  outward  into  an  unbounded  medium.  Therefore  in  summary,  the 
eigenfunctions  of  a  sphere  are  Pn( cos  d)  and  hXia\kr),  and  the  pressure  at  an 
arbitrary  point  may  be  wntten  in  terms  of  these  eigenfunctions  as 


P(r,0)  -  £  \Pn  (cost»h^]  (kr) 
n»o 
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